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Abstract 

In Gauss-Bonnet braneworld cosmology, the Friedmann equation of our four-dimensional 
universe on 3-brane is modified in a high energy regime (Gauss-Bonnet regime), while the 
standard expansion law is reproduced in low energies (standard regime). We investigate the 
Gauss-Bonnet braneworld cosmological effect on the thermal relic density of cold dark matter 
when the freeze-out of the dark matter occurs in the Gaugss-Bonnet regime. We find that the 
resultant relic density is considerably reduced when the transition temperature, which connects 
the Gauss-Bonnet regime with the standard regime, is low enough. This result is in sharp 
contrast with the result previously obtained in the Randall-Sundrum braneworld cosmology, 
where the relic density is enhanced. 



1 Introduction 



Recent various cosmological observations, in particular, the Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite p], have established the ACDM cosmological model with a great accu- 
racy and the relic abundance of the cold dark matter is estimated as (68 % CL. uncertainties) 

Vt CDM h 2 = 0.1131 ± 0.0034. (1) 

However, to clarify the identity of the dark matter particle is still a prime open problem in 
particle physics and cosmology. Since the Standard Model (SM) has no suitable candidate for 
the cold dark matter, the observation of the dark matter suggests new physics beyond the SM 
in which the dark matter particle is naturally provided. Many candidates for dark matter have 
been proposed in various new physics models, for example, the neutralino in supersymmetric 
models is one of the promising candidates [2]. 

Among several possibilities, the dark matter as the thermal relic may be the most plausible 
scenario since in this case, the relic abundance of the dark matter is insensitive to the history 
of the early universe before the freeze-out time of the dark matter, such as the mechanism of 
reheating after inflation etc. This scenario allows us to estimate the dark matter relic density 
by solving the Boltzmann equation [3], 

where Y = n/s is the yield defined by the ratio of the dark matter number density (n) to 
the entropy density of the universe s = 0A39g if m 3 /x 3 , g* ~ 100 is the effective total number 
of relativistic degrees of freedom, and x = m/T is the ratio between the dark matter mass 
(m) and the temperature of the universe (T). The yield in equilibrium Yeq is written as 
Y EQ = O.U5{g/g^ 2 e~ x for x > 3 with g being the degrees of freedom of the dark matter. In 
the standard cosmology, the Hubble parameter is described by the total energy density of the 
universe (p = (tt 2 /30)g*(m/x) 4 ) through the Friedmann equation of the form 

where Mpi = 1.22 x 10 19 GeV is the Planck mass. Using explicit formulas presented above, the 
Boltzmann equation can be rewritten into the form 

with an ^-independent constant A = xs/H = 0.26^/g^mMpi. 
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Solving the Boltzmann equation, we can obtain the thermal relic abundance of the dark 
matter at the present universe. As is well known, an approximate formula of the solution to 
the Boltzmann equation (for the S-wave annihilation process) can be described as 



where Xd = m/Td with the decoupling temperature Tj. It is useful to express the relic density 
in terms of the ratio of the dark matter density to the critical density, Qh 2 = msoY(oo)h 2 / p c , 
where p c = 1.1 x l(T 5 /i 2 cm" 3 , h = 0.7lj£{£ and s = 2900 cm" 3 [3]: 



In a typical dark matter scenario such as the WIMP (weakly interacting massive particle) 
scenario, Xd ~ 23 and thus, the dark matter relic density is controlled only by its annihilation 
cross section (av). 

Recently, the braneworld models have attracted lots of attention as a novel higher dimen- 
sional theory. In these models, it is assumed that the standard model particles are confined on 
a "3-brane" while gravity resides in the whole higher dimensional spacetime. The braneworld 
cosmology based on the model first proposed by Randall and Sundrum (RS) [3], the so-called 
RS II model, has been intensively investigated [5]. It has been found that [6] the Friedmann 
equation in the RS cosmology leads to a non-standard expansion low in high energies, while 
the standard cosmology is reproduced in low energies. 

Since the Hubble parameter is involved in the Boltzmann equation, the dark matter relic 
abundance depends on the expansion law of the early universe at freeze-out time. Therefore, if 
our universe is higher dimensional and obeys the non-standard expansion law at the freeze-out 
time, the resultant relic abundance of the dark matter will be altered. The RS braneworld 
cosmological effect on the dark matter physics has been investigated in detail [7J [8j [9j [10], in 
particular, it has been shown that the resultant dark matter relic abundance is considerably 
enhanced. In the same context of the RS braneworld cosmology, other cosmological issues such 
as leptogenesis [11] and the cosmological gravitino problem in super symmetric models [12] have 
also been examined. 

In this letter, we investigate the non-standard cosmological effect on the dark matter relic 
abundance in the context of the Gauss-Bonnet (GB) braneworld cosmology. Once the Gauss- 
Bonnet term is added in the RS braneworld model, the Friedmann equation in high energies is 
modified from the one in the RS braneworld cosmology. We will show that the GB braneworld 
cosmological effect works to reduce the thermal relic abundance of the dark matter. This is in 
sharp contrast with the RS braneworld cosmological effect. 




(5) 



tth 2 



1.07 x 10 9 x d GeV 



(6) 
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2 Relic density under the modified Friedmann equation 



Let us begin with parameterizing a modified Friedmann equation in a certain class of braneworld 
cosmological models such as 

H = H st F(x t /x), (7) 

where H st is the Hubble parameter in the standard cosmology, and the function F denotes 
the modification of the Friedmann equation in braneworld models. Phenomenologically viable 
models should reproduce the standard cosmology in low energies, so that we impose a condition: 
F(x t /x) = 1 for x t jx < 1. Here, x t = T t /m and we call T t "transition temperature" at which 
the modified expansion law shifts into the standard cosmological one. The model-independent 
constraints on T t is given by the success of the Bing Bang Nucleosynthesis (BBN) in the 
standard cosmology and the transition should complete before the BBN era, T t > 1 MeV0- For 
concreteness, we assume 

F(x t /x)={^)\ (8) 

for Xtjx > 1 with a constant 7. We can see that this parameterization is in fact a good 
approximation for Friedmann equations obtained in known braneworld cosmological models. 
For example, 7 = 2 corresponds to the RS braneworld cosmology, while 7 = —2/3 to the GB 
braneworld cosmology as we will see later. Here we keep 7 as a free parameter to make our 
discussion applicable to general braneworld models (if such models exist). 

First we approximately solve the Boltzmann equation with the modified Friedmann equation 
in the non-standard cosmology. Eq. (j4]) is modified into 



dY X ( (av) 



(Y 2 ~Y'), (9) 



dx x 2 \F(xt/x) J E ® ' 

from which we can understand that the effect of the modified Friedmann equation is equivalent 
to modify the annihilation cross section in the standard cosmology; 

M-C-J=y = <™>C-Y a") 



F(x t /x)J \x t 

in the era x t /x > 1 which we are interested in. Eq. implies that the braneworld cosmological 
effect enhances (reduces) the thermal relic abundance of dark matter for 7 > (7 < 0). 



1 The precision measurements of the gravitational law in sub-millimeter range lead to much more stringent 
constraint, for example, Tt > 1 TeV for the RS braneworld model [Jj. However, this constraint is, in general, 
quite model-dependent and can be moderated in some extended models [13] ■ In this paper, we consider only 
the model-independent BBN constraint on the transition temperature. 
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For simplicity, we parameterize the thermal average of the annihilation cross section times 
the relative velocity as (av) = a n x~ n with a (mass dimension 2) constant o n and an integer n 
(n = and 1 correspond to S-wave and P-wave processes, respectively). At the early time, the 
dark matter particle is in the thermal equilibrium and Y tracks Yeq closely. To begin, consider 
the small deviation from the thermal distribution A = Y — Yeq <C Yeq. The Boltzmann 
equation leads to 



x 2 {x t /x) 1 - 



'EQ 



A ~ - )—L—L dx ~ Xt x 2+ "~ 7 (11) 
\a n x- n (2Y EQ + A) 2\a n ' v ; 

where we have used an approximation formula cIYeq/cIx ~ —Yeq. As the temperature decreases 
or equivalently x becomes large, the deviation grows since Yeq is exponentially dumping. Even- 
tually the decoupling occurs at Xd roughly evaluated as A(xd) — Y(xd) — YEQ^Xd). At further 
low temperature, A ~ Y ^> Yeq is satisfied and Fjg term in the Boltzmann equation can be 
neglected so that 

^ = 2 A 2 , (12) 

dx x t 

and the solution is formally given by 

A{x) = AjxJ) + ( 7 -n-l)x 7 ^ n 1 " ^ U ^ ' (13) 

For the S-wave process (n = 0), for example, the well-known result in the standard cos- 
mology (corresponding to 7 = 0), l/F(oo) ~ Xcr /xd, is obtained. When we take 7 = 2, our 
analysis here is the same as the one in [7j for the RS braneworld cosmology. For n = 0, for 
example, A(x) _1 is continuously growing and this growth stops at X — Xt 5 SO that the resultant 
relic abundance has been found to be l/y(oo) ~ \(jQX t . Then, we obtain the ratio of the energy 
density of the dark matter in the RS braneworld cosmology (ft(Rs)) to the one in the standard 
cosmology (0( s )) such that 

%a * (*_) , (i4) 

"(s) \Xd{s) J 

where xm s ) is the decoupling temperature in the standard cosmology. The relic abundance is 
enhanced by the RS braneworld effect for x t > Xd( s ), while it should saturate to the standard 
result for x t < Xd( s )- 

On the other hand, if 7 < 0, we arrived at the result which is in sharp contrast with the 
one in the RS braneworld cosmology. For simplicity, we consider only n = in the following 
and the analysis for n > is straightforward. In this case, we obtain 

x d , (15) 



F(oo) I-7 \x d 
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and the ratio of the energy density of the dark matter in this braneworld cosmology (fta)) to 
the one in the standard cosmology (f2( s )) is evaluated as 



(16) 

'-{a) \ x t J \Xd(s)J 

Thus, the resultant relic density is reduced in the case for 7 < and Xd < x t . This is nothing 
but the case that happens in the GB braneworld cosmology. 

3 Gauss-Bonnet braneworld cosmological effect 

Motivated by string theory considerations, it is a natural extention to add higher curvature 
terms to the bulk gravity action in the RS braneworld model [H]. Among possible terms, the 
Gauss-Bonnet invariant is of particular interests in five dimensions, since it is a unique nonlinear 
term in curvature which yields second order gravitational field equations. The five-dimensional 
gravitational action with the GB invariant is given by 

S = lh I d 5 x^gT [-2A 5 + K + a (TZ 2 - A1Z ab 1l ab + TZ abcd n abcd )} 
Zni J 



d X V /3 #4 (CT + C-matter) , (17) 

' brane 

where k\ = 87r/M| with the five-dimensional Planck scale M 5 , a > is the brane tension, and 
A 5 < is the bulk cosmological constant. The RS model is recovered in the limit a — ► 0. 

Imposing a Z 2 parity across the brane in an anti-de Sitter bulk and modeling the matters 
on the brane as a perfect fluid, the modified Friedmann equation on the spatially flat brane has 
been found as 



4(p + ml) = 2 H j 1 + ^ (3 - + 2/3^ ) , (18) 



where (3 = Aafi 2 = 1 — a/1 + 4aA 5 /3 and m a = cr 1 / 4 . There are four free parameters, k 5 , m a , 
li and P, corresponding to the original free parameters, k$, a, A 5 and a, which are constrained 
by phenomenological requirements. To reproduce the Friedmann equation of the standard 
cosmology with zero cosmological constant in the limit H 2 /fi 2 < 1, we find two relations 
among the parameters: 

K *=wr it^ 5 - (19) 
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The modified Friedmann equation can be rewritten in the useful form [T7] 

'2* 



tt2 V 2 
H = 1 



;i-/3)cosh(^) -1 



p + m A a = sinh x, (20) 
where x is a dimensionless measure of the energy density and 



( «-f_ y (i - ay 
m "-''— *4 2 <fWTW- (21) 

Here we have used Eq. f fT9|) to eliminate k 5 in the last equality In the same way, we express 
m„ as 



'k 2 4 \i + p 



K = 2 ^ fr4 • (22) 



The evolution of the GB braneworld cosmology is characterized by the two mass scales, m a 
and ma. Expanding Eq. (120]) with respect to x> we find three regimes for m a > rn a : The GB 
regime for p ^> m^, 

H 2 * (^-»4pY 3 , (23) 



4/? 



the RS regime for ^> p ^> in 4 



H 2 * ^p\ (24) 



and the standard regime for m A a ^> p, 



# 2 - fp- (25) 



Since we are interested in the GB regime, let us simplify the evolution of the universe by 
imposing the condition m a = m a , which leads to 

3/5 3 - 12p 2 + 15/5 - 2 = (26) 

and hence, f3 = 0.151. In this case, the RS regime is collapsed and there are only two regimes in 
the evolution of the universe. Applying the parameterization to the non-standard Friedmann 
equation in Eqs. ((Tj) and (JHJ), 

H = H st (*) 16 = H st (-) 2 3 (27) 



P J \ X tJ 
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Figure 1: The Hubble parameter as a function of p in the GB braneworld cosmology. The solid 
line shows the exact formula of the Friedmann equation, while the dashed (in red) corresponds 
to our approximation formula. Here, we have taken a unit — 



Pt = V-A-^ i ) 3^ 3 '«, (2c 



for p > p t or equivalently x < x t , while H = H st for p < p t , where 

27 A + ZAV^,,^ 

16 v P J 4 

In Fig. 1, we show that our approximation for the Friedmann equation is in fact a good ap- 
proximation to the exact form in Eq. ( 1201) . 

Now we are ready to see the GB braneworld cosmological effect on the dark matter relic 
abundance. Eq. (|27|) means that the modified Friedmann equation in the GB braneworld 
cosmology corresponds to 7 = —2/3, and thus we find 



Q (gb) _ 5 fxA 2/i 



x d 



(29) 



^(s) 3 \X d J \Xd{s), 

from Eq. (|T6l) . Therefore, the GB braneworld cosmological effect reduces the relic density and 
the reduction rate is controlled by the transition temperature. This is our main result. For the 
WIMP dark matter, the typical value of the decoupling temperature is Xd ~ 23 and this is not 
so much changing even under the non-standard Friedmann equation. Thus, we expect 

Oca*) „ 5 _ (30) 



f2( s ) 3 V23/ 

For example, fi(GB)/^( s ) — 0.25 for x t = 400. 

Finally, let us check that our analytic result given above is a good approximation for the 
results from the numerical solution of the Boltzmann equation. Fixing the dark matter mass 
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Figure 2: Numerical solutions of the Boltzmann equation in the Gauss-Bonnet braneworld 
cosmology (solid line) and in the standard cosmology (dashed line). The dotted line corresponds 
to Yeq. The GB braneworld cosmological effect reduces the resultant relic abundance from the 
one in the standard cosmology. 

m = 100 GeV, (av) = 10~ 5 /m 2 and Xt = 400, we numerically solve the Boltzmann equation 
with the Friedman equations in the standard cosmology and the GB braneworld cosmology, 
respectively. The numerical results are depicted in Fig. 2. Here, we obtain 0(gb)/^(«) — 0.26, 
which is very close to our previous result from analytic formulas. 

4 Conclusions and discussions 

We have investigated the thermal relic density of the cold dark matter in the braneworld cosmol- 
ogy, in particular, the Gauss-Bonnet braneworld cosmology, which is a natural extension of the 
Randall- Sundrum braneworld model to include the higher curvature terms. We have modeled 
the modified Friedmann equation in such a way applicable to general braneworld cosmological 
models and analytically solved the Boltzmann equation under some approximation. Applying 
this result to the Gauss-Bonnet braneworld cosmology, we have found that the resultant relic 
density of the dark matter is considerably reduced from the one in the standard cosmology, when 
the freeze-out occurs well in the Gauss-Bonnet regime in the evolution of the universe. This 
conclusion is in sharp contrast with the result in the Randall-Sundrum braneworld cosmology, 
where the relic density is enhanced by the braneworld cosmological effect. 

It is worth applying our results in this paper to concrete dark matter models which have 
been investigated in the standard cosmology. For supersymmetric models with the neutralino 
dark matter, the RS braneworld cosmological effect was analyzed [H]. It has been shown that 
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the allowed parameter region for the neutralino dark matter consistent with the observed dark 
matter density is dramatically modified from the one in the standard cosmology and eventually 
disappears as the transition temperature is lowered. This is because the RS braneworld cosmo- 
logical effect enhances the dark matter relic density while supersymmetric models, in particular, 
the constrained minimal supersymmetric SM, tend to predict an over-abundance of neutralino 
dark matter. On the other hand, the Gauss-Bonnet braneworld cosmological effect reduces 
the relic density of the dark matter and therefore can enlarge the allowed parameter region 
in supersymmetric models. A similar effect has been discussed in the scalar-tensor cosmology 
[18j. This enlargement of the cosmologically allowed parameter region has an impact on the 
sparticle search at the Large Hadron Collider (LHC). As have been investigated in [TTj IT2] for 
the RS braneworld cosmology, it would be interesting to consider other cosmological issues also 
in the GB braneworld cosmology. We leave these subjects in future works. 
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